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Abstract 

The parameters Xi.2> which measure nonfactorizable soft gluon contributions to hadronic weak 
decays of mesons, are updated by extracting them from the data of D, B — > PP, VP decays (P: 
pseudoscalar meson, V: vector meson). It is found that \2 ranges from —0.36 to —0.60 in the 
decays from D — > Kit to D + — ► <j)7r + , D — > K*tt, while it is of order 10% with a positive sign in 
B — > ipK, Dtt, D*tt, Dp decays. Therefore, the effective parameter ai is process dependent in 
charm decay, whereas it stays fairly stable in B decay. This implies the picture that nonfactorizable 
effects become stronger when the decay particles become less energetic after hadronization. As for 
D, B — > VV decays, the presence of nonfactorizable terms in general prevents a possible definition of 
effective a\ and ai- This is reenforced by the observation of a large longitudinal polarization fraction 
in B — > ipK* decay, implying 5-wave dominated nonfactorizable effects. The nonfactorizable term 
dominated by the S- wave is also essential for understanding the decay rate of B~ — > D*°p~ . It is 
found that all nonfactorizable effects /Af K * , A^/Af p , A\^ /Af D * (nf standing for nonfactor- 
ization) are positive and of order 10%, in accordance with Xi(B — ► D(D*)ir(p)) and X2(B — ► ipK). 
However, we show that in D — > K* p decay nonfactorizable effects cannot be dominated by the S- 
wave. A polarization measurement in the color- and Cabibbo-suppressed decay mode D + — > (pp + 
is strongly urged in order to test if plays a more pivotal role than A± /A\ in charm decay. 
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1. Introduction 



It is customary to assume that two-body nonleptonic weak decays of heavy mesons are 
dominated by factorizable contributions. Under this assumption, the spectator meson decay 
amplitude is the product of the universal parameter at (for external ^-emission) or a 2 (for 
internal jy-emission) , which is channel independent in D or B decays, and hadronic matrix 
elements which can be factorized as the product of two independent hadronic currents. The 
universal parameters a\ and a 2 are related to the Wilson coefficient functions c\ and C2 by 

1 1 

a\ = ci + — c 2 , a 2 =c 2 + — ci, (1) 

with iV c being the number of colors. It is known that the bulk of exclusive nonleptonic charm 
decay data cannot be explained by this factorization approach [1]. For example, the predicted 
ratio of the color-suppressed mode D° — > K°7t and color-favored decay D° — ► K~n + is in 
violent disagreement with experiment. This signals the importance of the nonfactorizable 
effects. 

The leading nonfactorizable contribution arises from the soft gluon exchange between 
two color-octet currents 

O c = ^lA^feA^), (2) 

where (qi\ a q2) stands for gi7 M (l — 7s)A a g2- For M — > PP, VP decays (P: pseudoscalar me- 
son, V: vector meson), the nonfactorizable effect amounts to a redefinition of the parameters 
ai and a 2 [2], Q 

ax -> ci + c a(^ +Xi), a 2 -> c 2 + ci(^- + X2), (3) 

where Xi an d X2 denote the contributions of O c to color-favored and color-suppressed decay 
amplitudes respectively relative to the factorizable ones. For example, for Df — > 07r + , D + — > 
</)7r + decays, 



Xl (£>+-^r + ) 



+ \ _ W" 12 



7r+|i(uA a d)(sA a c)|D+) 



(^r+|(tid)(ffc)|D+> 



m + . +x (07r+|i(^c)( g -A a g )|D+) 

(07T + 1 («C) (SS)P + )/ 

The subscript / in Eq.(4) denotes a factorizable contribution: 

(^ + \(ud)(sc)\D+) f = 2mMe*-p Ds )A^(ml), 
(<P7T + \(uc)(ss)\D+) f = m<P U{e*.p D )F^{w?X (5) 



1 Note that our definition of \i an d X2 is different from r\ and T2 defined in [3] by a factor of 2. 
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where is the polarization vector of the meson, and we have followed Ref.[4] for the 
definition of form factors. The nonfactorizable contributions have the expressions 

(4>it + \\{u\ a d){s\ a c)\Dt) = 2mMe*-p Ds )A^(ml), 
(<P^(u\ a c)(s\ a s)\D+) = m 4 U(e*-p D )F? f (ml), (6) 
with the superscript nf referring to nonfactorizable contributions. It is clear that 

That is, x simply measures the fraction of nonfactorizable contributions to the form factor 
under consideration. 

Although we do not know how to calculate xi and Xi from first principles, we do anticipate 
that [3] 

\ X (B - PP)\ < \ X (D - PP)\ < \ X (D - VP)\, (8) 

based on the reason that nonperturbative soft gluon effects become more important when the 
final-state particles move slower, allowing more time for significant final-state interactions 
after hadronization. As a consequence, it is obvious that ai :2 are in general not universal and 
that the rule of discarding 1/N C terms [5], which works empirically well in D — > Kit decay, 
cannot be safely extrapolated to B — > Dir decay as \x(B — > Dtt)\ is expected to be much 
smaller than \x(D — > Kir) ~ — 1| (the cm. momentum in D — > Kir being 861 MeV, to be 
compared with 2307 MeV in B — > D7r) and hence a large cancellation between 1/N C and 
x(.B — > £)7r) is not expected to happen. The recent CLEO observation [6] that the rule of 
discarding 1/N C terms is not operative in B — > D(D*)ir(p) decays is therefore not stunning. 
Only the fact that x(B — > -D7r) is positive turns out to be striking. 

Unlike the PP or VP case, it is not pertinent to define Xi,2 for M — > VV decay as its 
general amplitude consists of three independent Lorentz scalars: 

A[M(p) - V 1 (e 1 , Pl )V 2 (e 2 ,p 2 )} oc ^(A 1 )<(A 2 )(i 1 ^ + A 2 p»p» + ^e^W 2/3 ), (9) 

where A 1 , A 2 , V are related to the form factors A 1 , A 2 and V respectively. Since a priori 
there is no reason to expect that nonfactorizable terms weight in the same way to S-, P- 
and D-waves, namely A™ f /A 1 = A 2 lf /A 2 = V nf /V, we thus cannot define Xi an d Xi- Conse- 
quently, it is in general not possible to define an effective a x or a 2 for M — > VV decays once 
nonfactorizable effects are taken into account [7] . In the factorization approach, the fraction 
of polarization, say T L /T (L: longitudinal polarization) in B — > ipK* decay, is independent 
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of the parameter a 2 . As a result, if an effective a 2 can be defined for B — > ipK*, it will lead to 
the conclusion that nonfactorizable terms cannot affect the factorization prediction of Yl/Y 
at all. It was realized recently that all the known models in the literature in conjunction 
with the factorization hypothesis fail to reproduce the data of Yl/Y or the production ratio 
Y(B — > ipK*)/Y(B — > ipK) or both [8,9]. Evidently, if we wish to utilize nonfactorizable ef- 
fects to resolve the puzzle with Yl/Y, a key ingredient will be the nonexistence of an effective 
a 2 for B -> ipK*. 

In short, there are two places where the factorization hypothesis can be unambiguiously 
tested: (i) To extract the parameters ai and a 2 from the experimental measurements of 
M — ^ PP, VP to see if they are process independent, (ii) To measure the fraction of 
longitudinal polarization in M — > VV decay and compare with the factorization prediction. 
Any failure of them will indicate a breakdown of factorization. 

The purpose of the present paper is threefold, (i) The parameters xi an d X2 have been 
extracted in Ref.[3] (see also [10]). Here we wish to update the values of Xi,2 using the q 2 
dependence of form factors suggested by QCD-sum-rule calculations and other theoretical 
arguments, (ii) It was recently advocated by Kamal and Sandra [7] that the assumption that 
in B — > ipK* decay the nonfactorizable amplitude contributes only to S-wave final states, 
namely 7^ 0, A 2 ^ = V n ^ = 0, will lead to a satisfactory explanation of the data of 
Y(B — > vpK*)/Y(B — > ipK) and Yl/Y. We would like to show that this very assumption is 
also essential for understanding the ratio B(B~ — > D*°p~)/B(B° — > D* + p~), which cannot 
be explained satisfactorily in previous work, (iii) Contrary to the B meson case, we will 
demonstrate that the assumption of S'-wave dominated nonfactorizable terms does not work 
in D — > VV decay. 

2. Nonfactorizable contributions in D, B — > PP, VP decays 

Because of the presence of final-state interactions (FSI) and the nonspectator contribu- 
tions (W^-exchange and W^-annihilation), it is generally not possible to extract the nonfac- 
torization parameters xi,2 except for a very few channels. Though color-suppressed decays, 
for example, D° — > K° \K*®)n G '(p°) are conventionally classified as Class II modes [11], color- 
flavored decay D° — > K~tt + will bring some important contribution to D° — > K°tc via 
FSI. This together with the small but not negligible W -exchange amplitude renders the de- 
termination of 02 from D° — > K°tc° impossible. Therefore, in order to determine a± and 
especially a 2 we should focus on the exotic channels e.g. D + — > K°tt + , tt + 7t , and the decay 
modes with one single isospin component, e.g. D + — > 7r + 0, — > 7r + 0, where nonspectator 
contributions are absent and FSI are presumably negligible. 
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We next write down the relations between xi,2 and form factors 



F DK (ml)' 7 F^\ 



m 



K) 



Xl (D - K*n) = ^fl , X 2(D - K*n) ^ ^ 



^r p)= *&>, MD ^ p) .gm, (l0 , 



It is clear that only the three form factors F , F 1 and A entering into the decay amplitudes 
of M — > PP, VP. A consideration of the heavy quark limit behavior of the form factors 
suggests that the q 2 dependence of Pi (A 2 ) is different from that of P (A and A{) by an 
additional pole factor [12]. Indeed, QCD-sum-rule calculations have implied a monopole 
behavior for Fi(q 2 ) [13-16] and an approximately constant F [15]. With a dipole form 
factor A 2 , as shown by a recent QCD-sum-rule analysis [16], we will thus assume a monopole 
behavior for Aq. 

Unlike the decays D + — > 7r + 0, P+ — > 7r + which are described by a single quark diagram, 
we cannot extract xi,2 from the data of D + — > PJ°7r + , K°p + , K*°tc + alone without providing 
further information. For example, the decay amplitude of D + — > PJ°7r + reads 

A(P+ - P;°7r + ) = ^V: s V ud [a,{ml - m 2 K )f n F DK (ml) + a 2 (m 2 D - ml)f K F ^(m 2 K )}, (11) 

which consists of external lU-emission and internal lU-emission amplitudes. We will therefore 
make a plausible assumption that Xi ~ X2 so that x(F> — > Kn) can be determined from 
the measured rate of D + — > K°ir + . Since the extraction procedure is already elucidated in 
Ref.[3], here we will simply present the results (only the central values being quoted) followed 
by several remarks 

X2 (D -> Kn) ~ -0.36, 

X 2 (D ^ K*n) ~ -0.61, (12) 
X 2 (P + ^0vr + ) ~ -0.44, 

where we have used the following quantities: 

ci(m c ) = 1.26, c 2 (m c ) = -0.51, 

f n = 132 MeV, f K = 160 MeV, f K , = 220 MeV, = 237 MeV, 

F° K {0) = F? K (0) = 0.77 ± 0.04 [17], P D "(0) = P 1 D7r (0) = 0.83 [18], (13) 

A? K *(0) = 0.61 ±0.05, A^ x *(0) = 0.45 ±0.09 [17], < x *(0) = 0.70, 
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and the Particle Data Group [19] for the decay rates of various decay modes. 

Several remarks are in order, (i) As pointed out by Soares [10], the solutions for x are n °t 
uniquely determined. For example, the other possible solution for XiiP — * Kn) is —1.18. 
To remove the ambiguities, we have assumed that nonfactorizable corrections are small 
compared to the factorizable ones, (ii) Assuming Aq P (0) = Aq K *(0), we find from the decay 
D + — > K°p + that x(D — > Kp) rs —1.5, wich is unreasonably too large. We do not know how 
to resolve this problem except for noting that thus far there is only one measurement of this 
decay mode [20]. (iii) To determine Xi{P ) t ~^ 0tt + ) requires a better knowledge of the form 
factor Aq s ^ and the branching ratio of — > (f)7i + . Unfortunatly, a direct measurement of 
them is still not available. Assuming A%'*(0) « A° K * (0) and B(D+ -> 0tt+) = (3.5 ±0.4)% 
[19], we get xi(Df — > 07r + ) —0.60. So in general nonfactorizable terms are process or 
class dependent, and satisfy the relation \x(D — > PP)\ < \x{D ~^ VP)\ as expected, (iv) 
Since X2(D — * Kit) is close to — |, it is evident that a large cancellation between 1/N C and 
X2(D — > Kir) occurs. This is the dynamic reason why the large- A^, approach operates well 
for D — > Ktt decay. However, this is no longer the case for D — > VP decays. The predicted 
branching ratios in 1/N C expansion are 

B{D + -> K*\ + ) = 0.3%, B{D + -> #V) = 16%, 

£(£+ -> = 17%, B{D + -> 07T+) = 0.4%, (14) 

to be compared with data [19] 

^ ^*°7r+) cxpt = (2.2 ± 0.4)%, £(D + -> X p + )expt = (6.6 ± 2.5)%, 

£ (jD + ^ K*V + )ex P t = (4.8 ± 1.8)%, £(£+ -> 07r+) cxpt = (0.67 ± 0.08)%. (15) 

Consider the decay D + — > K*°n + as an example. Its amplitude is given by 

A (D + ^ R*\ + ) = V2G F V c * s V ud [aiUm K *A° K *(ml) + a 2 f K ,m K ,F^{m 2 K ,)}. (16) 

Since the interference is destructive and fK*Ff )7r > !-kA® k * , a large \a 2 \ is needed in order to 
enhance the branching ratio of D + — > K*°ir + from 0.3% to 2.2%. (Note that a\ is relatively 
insensitive to the nonfactorizable effects.) This in turn implies a negative (jj- + X2) and 
hence Xi{D — > K*n) < — |. Therefore, we are led to conclude that the leading 1/N C 
expansion cannot be a universal approach for the nonleptonic weak decays of the meson. 
However, the fact that substantial nonfactorizable effects which contribute destructively 
with the subleading 1/N C factorizable contributions are required to accommodate the data 
of charm decay means that, as far as charm decays are concerned, the large- N c approach 
greatly improves the naive factorization method in which xi,2 = 0; the former approach 
amounts to having a universal nonfactorizable term xi,2 = — 1/N C . 
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We next turn to B — > D(D*)ir(p) decays. Though both nonspectator and FSI effects 
are known to be important in charm decays, it is generally believed that they do not play a 
significant role in bottom decays as the decay particles are moving fast, not allowing adequate 
time for FSI. This gives the enormous advantage that it is conceivable to determine a\ and 
a 2 separately from B — > D(D*)ir(p) decays. Using the heavy-flavor-symmetry approach for 
heavy-light form factors and assuming a monopole extrapolation for Fx, A , Ai, a dipole 
behavior for A 2 , V, and an approximately constant F , as suggested by QCD-sum-rule 
calculations and some theoretical arguments [21], we found from the CLEO data that [21] 



a 2 (B - DWtt(p)) 



1.01 ±0.06, 
0.23 ±0.06. 



;i7) 



Taking Ci(m&) = 1.11 and c 2 (rrib) = —0.26 leads to 

Xi(B -> D { *\(p)) ~ 0.05 , xi[B -> D w 7r(p)) ~ 0.11 



18) 



Since (^- + ^1,2) = (01,2 — ci,2)/c 2 ,i and | C2 1 << |ci|, it is clear that the determination 
of Xi is f ar more uncertain than X2- it is very sensitive to the values of a±, c\ and c 2 . 
We see from (18) that nonfactorizable effects become less important in B decays, as what 
expected [see (8)]. However, a positive X2(B — > D(D*)ir(p)), which is necessary to explain 
the constructive interference in B~ — > D° (D* )tt~ (p - ) decays, appears to be rather striking. 
A recent light cone QCD-sum-rule calculation [22] following the framework outlined in [23] 
fails to reproduce a positive X2(B — > Dtt). This tantalizing issue should be resolved in the 
near future. 



For B — > if)K decays, we found [21] 



a 2 {B~ -> ipK~ 



0.235 ±0.018 



a 2 (B° -> ipK c 



0.192 ±0.032. 



(19) 



The combined value is 



a 2 (B -> ipK) = 0.225 ± 0.016 , (20) 

where its sign should be positive, as we have argued in [21]. (It was advocated by Soares 
[10] that an analysis of the contribution of B — > ipK to the decay B — > Ki + i~ can be used 
to remove the sign ambiguity of a 2 .) It follows that 

F nf (m 2 ) 

X*{B^1>K)= l K \J -0.10, (21) 

2 Contrary to the charmed meson case, the variation of 0^2 from B — > Dtt to D*ir and Dp decays is 
negligible (see Table IV of [21]). 
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which is in accordance with x%{B — » 0"Jr()5)). 

Finally, it is very interesting to note that, in contrast to charm decays, the large- iV c 
approach is even worse than the naive factorization method in describing B — > D(D*)7r(p) 
decays as X2(B — > D(*'ir(p)) is small but positive. 

3. Nonfactorizable contributions in B — > ipK*, D*p decays 

As stressed in the Introduction, in general one cannot define Xi,2 and hence an effective 
dip for M — > VV decays unless the nonfactorizable terms weight in the same manner 
in all three partial waves. It was pointed out recently that there are two experimental 
data, namely the production ratio R = Y(B — > ipK*)/Y(B — > ipK) and the fraction of 
longitudinal polarization Yl/Y in B — > ipK*, which cannot be accounted for simultaneously 
by all commonly used models within the framework of factorization [8,9]. The experimental 
results are 

R= 1.74 ±0.39 [6], ^=0.78 ±0.07, (22) 
where the latter is the combined average of the three measurements: 



r 



0.97 ± 0.16 ± 0.15, ARGUS [24]; 

0.80 ±0.08 ±0.05, CLEO [6]; (23) 
U ~' " s 1 0.66 ±0.1018$, CDF [25]. 



Irrespective of the production ratio R, all the existing models fail to produce a large 
longitudinal polarization fraction [8,9]. This strongly implies that the puzzle with Yl/Y can 
only be resolved by appealing to nonfactorizable effects. However, if the relation A" J ' / ' A\ = 
A n 2 s jA 2 = V nf /V holds, then an effective ci2 can be defined for B — > ipK* and the prediction 
of Yl/Y will be the same as that in the factorization approach as the polarization fraction is 
independent of a 2 . Consequently, nonfactorizable terms should contribute differently to S-, 
P- and D-w&ve amplitudes if we wish to explain the observed Yl/Y. 

The large longitudinal polarization fraction observed by ARGUS and CLEO suggests 
that the decay B — > ipK* is almost all S-w&ve. To see this, we write down the B — > ipK* 



3 An interesting observation was made recently in [26] that the factorization assumption in B — » t/jK(K*) 
is not ruled out and the data can be accommodated by the heavy-flavor-symmctry approach for heavy-light 
form factors provided that the Ai(q 2 ) form factor is frankly decreasing. To our knowledge, a decreasing 
A\ with q 2 is ruled out by several recent QCD-sum-rule analyses (see e.g. [16]). Using the same approach 
for heavy-light form factors but the q 2 dependence of form factors given in [21], we found that R = 1.84 
and Tl/T = 0.56 [21]. Evidently, the factorization approach is still difficult to explain the observed large 
polarization fraction. 
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amplitude 
A[B(p) 

with 



T2 



Cl 



V: s V ud (c2 + ^) Um^)el{K*) [A l9 ^ + A^f 



+ iVe^p la p 2p ], 



(24) 



Ai = (m B + m K *)A 1 



BK* , 



Ao = 



V 



(m B + m K *) 
2 

(m B + m K *) 



A 



BK* 



V BK \mi) 



1 + K 
1 + K 



A n 2 f H) 



v nf H) 

V BK \m 2 ) 



(25) 



and k = Ci/ (c 2 + §ci). It is easily seen that we will have an effective a 2 = c 2 + ci(| + X2) if 
the nonfactorizable terms happen to satisfy the relation A™* /Ai = A 2 j A 2 = V n f /V = Xi- 
The decay rate of this mode is of the form 



T(B -> 4>K*) oc (a - bx) 2 + 2(1 + c 2 y 2 ), 



(26) 



where 



m B — — mr K * 



x 



^(^)M B/ K) 



b = 



y 



2m 2 B p 



2m B p c 



m^m K *(m B + m K *) 2 ' 
V BK \ml) + KV n f(m%) 



(m B + m K *) 2 ' 



(27) 



Af^(mJ) + «A?'(mJ)' 
with p c being the cm. momentum. The longitudinal polarization fraction is then given by 

(a — bx) 2 



r 



(a - bx) 2 + 2(1 + c 2 y 2 ) 



(28) 



If the decay is an almost S-wave, one will have T L /T ~ a 2 / (a 2 + 2) = 0.83 . Since k » 1, x 
(D-wave) and y (P-wave) can be suppressed by assuming that, as first postulated in [7], in 
B — > ipK* decay the nonfactorizable amplitude contributes only to S'-wave final states; that 
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IS. 



Af ^ 0, Al s = V nf = 0. 

The rational for this assumption is given in [7]. 

With the assumption (29), the branching ratio followed from (24) is 

Cl 



(29) 



B{B -> ijiK*) = 0.0288 



c 2 



[K-k) 2 + 2(^ + c y)] 



with 



2 \ ' 



V BK \m%) 



£ = i + «- 



where uses of \V r 



cb\ 



0.040 and r(S) = 1.52 x 10~ 12 s have been made. It follows that 
(a£ — bx) 2 



(30) 



(31) 



r 



(32) 



We use the measured branching ratio B[B — ► 
ratio (mfy/Af K *(mV), which is found to be 



<-fa;) 2 + 2(l + cV)' 

(0.172 ±0.030)% [6] to determine the 



Af K *(m 



0.08 



which we have used A\ 



BK* 



.m: 



0.41, AW(m 



(33) 



0.36, V BK \m%) = 0.72 [21] and 



discarded the other possible solution A\^ /Af K * = —0.22 for its "wrong" sign, recalling 
that Fi/Ff K is positive [cf. Eq.(21)]. The predicted longitudinal polarization fraction is 
Tl/T = 0.73, which is in accordance with experiment. 

The assumption of negligible nonfactorizable contributions to P- and D-w&ves also turns 
out to be essential for understanding the decay rate of B~ — > D*°p~ or the ratio R4 = 
B{B~ —> D*°p~)/B(B° —> D* + p~). The issue arises as follows. In Ref.[21] we have deter- 
mined a>i and a% from B — > Dn, D*ir, Dp decays and obtained a consistent ratio 02/01 from 



4 A different approach for nonfactorizable effects adopted in Ref.[27] amounts to A\ 
V nf ^ 0. It follows from Eq.(28) that 



nf - Af = and 



Tl 

r 



(a - bx) 2 



(a -bx) 2 + 2(1 + c 2 y 2 Y 



with y = (V BK {m 2 ^) + kV 1 ^ (m 2 b ))/Af K (rr&) and x being defined in (31). It is clear that in order to 
get a large longitudinal polarization fraction one needs a negative V n f /V ! Using the numerical values 
a = 3.164, b = 1.304, x = 0.89, we find (T L /T) max = 0.67. The prediction T L /T = 0.65 given by [27] is one 
standard deviation from experiment (22). 
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channel to channel: 0.24 ± 0.10, 0.24 ± 0.14, 0.21 ± 0.08 (see Table IV of [21]). Assuming 
factorization, we got a 2 ja\ — 0.34 ± 0.13 from B — > D*p decay, which deviates somewhat 
from above values. In the presence of S'-wave dominated nonfactorizable contributions, it 
is no longer possible to define an effective a\ and 02 for B — > D*p decay. Therefore, the 
quantities to be compared with are A™ jA\ in B — > D*p decay and % 2 in B — > Dn, D*n, Dp. 
A straightforward calculation yields (see [21] for the factorizable case) 

t(B~) ( Hi „H 2 

with 



* = imK 1+ *'ir + 'fir>> (34) 



H x = (ai-bx)(ai'-b'x') + 2(ii' + cc'yy'), 
H 2 = (ai'-b'x / ) 2 + 2(i /2 + c /2 y /2 ), 

m D *(m B + m p ) f D , Af p (m^.) c 2 + |ci 
m p (m B + m D *) f p Af D * (m 2 p ) c x + \c 2 ' 



e = 1 + 



c 2 Af(m 2 



ci + \c 2 Af D *[w 



2\ 



? = 1 + 



ci <(m 2 D , 



c 2 + | Cl Af'(rr£,)' 

where a, 6, c are obtained from a, b, c respectively in (27), x, y from x, y in (31) by 
replacing ip — > £>*, X* — > p, and b', c', x', y' are obtained from b, c, x, y respectively by 
replacing D* <-> p; for instance x' = Af p (m^,)Mf P ("ii?0- Assuming A? /A? D * ~ A? /Af p 
and fitting (34) to the experimental value i? 4 = (1.68 ± 0.35)% [6], we get 

<^~if^K)" a12 ^ (36) 

We see that the S-wave dominated nonfactorizable effect in B — > ipK* and B — > D*p decays 
is of order 10%, consistent with X2(B — > ipK) and %2(-6 — > D(D*)n(p)). 

4. Nonfactorizable contributions in D — > X*p decay 

We have shown in the previous section that S'-wave dominated nonfactorizable terms 
are needed to explain the large longitudinal polarization fraction observed in B — > 
and the ratio £>(£?~ — > D*°p~)/B(B° — > D* + p~). However, we shall see in this section 
that the assumption (29) is no longer applicable to D — > K*p decay. An experimental 
measurement of £>+ -> i?*°p + and L>° -> i?*°p° by Mark III [28] shows that (i) the decay 
D + — > K*°p + is a mixture of longitudinal and transverse polarization consistent with a pure 
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S-wave amplitude, f] and (ii) D° — > K*°p° is almost all transverse, requiring a cancellation 
between the longitudinal S-wave and Z)-wave. 



We first consider the decay D 
A(D + (p) 



K*°p + , whose amplitude is given by 



K*°(pi)p + ( P 2)) = ^V; s V ud e;(K*)e:(p)[A ig ^ + A^f + iVe^p XaPw ], (38) 



where 



A, 



ci + y ) f p m p (m D + m K *) ( 1 + 



c 2 



+ 



c 2 + — fK*mK*{m D +m p )[l + — -g- f 

3 / V c 2 + 3 c i Ai P (m 2 K 




(39) 



and A 2 (V) is obtained from Ax with the replacements Ai — > A 2 {A\ — > V), (tod + m^*) — > 
-2/(mc + mx*) and (m^ + m p ) — >■ -2/(mo + m p ). Since A^ /A^ K * and A^ /Af p are 
expected to be negative [see Eq.(12)], it is obvious that if nonfactorizable terms are dominated 
by the S'-wave, it will imply a more severe destructive interference in the S'-wave amplitude 
than in P- and D-w&ve amplitudes, in contradiction to the observation that this decay is 
almost all S-wave. The branching ratio is calculated to be 

2 



B(D + -> K*°p 



0.10 



ci 



(40) 



with the expressions of r/, if', ii{ 2 analogous to 77, if, iii >2 in (35). A fit of (40) to the 
Mark III data for the branching ratio (37) gives rise to (assuming A\ /Af K * ~ A^ /A^ p ) 



Af K *\ 



-0.98 



(41) 

which is uncomfortably too large. Q Moreover, the P-wave branching ratio is predicted to 
be 2.0 x 10 -2 , in disagreement with experiment [28] 



B{D + -> K*°p + )P-wave < 0.5 X 10' 

It thus appears to us that an almost S-wave D + — > K*°p + implies that 



A n 2 f 




ynf 


> 


A? 


A DK*(p) 


5 


yDK*{ P ) 


A DK*{p) 



(42) 



(43) 



5 The other measurement by E691 [29] disagrees severely with Mark III on the branching ratio 

f (4.8 ± 1.2 ± 1.4)%, Mark III [281; , , 

B(D+ ^ K* a p+) = I ) ' r / J ' 37 

v P 7 \ (2.3 ± 1.2 ±0.9)%, E691 [29]. v ; 

Recall that model calculations tend to give a very large branching ratio of 17% [see Eq.(14)]. 

6 A fit to the E691 measurement (37) for the branching ratio yields an even larger value: A™^ /A® K ~ 
A\ ! IA^ P ^ -1.41. 
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Taking A" f /A 1 = A\ s /A 2 = V nf /V = x(D -> K*p) as an illustration, we obtain 

X (D -> A» « -0.65 (44) 

and B(D + — > -^*°p + )p- W avc = 2.0 x 10~ 3 , which are certainly more plausible than before. 

Another indication for the failure of the S-wave dominated hypothesis for nonfactorizable 
effects comes from the decay D° — > K*°p°, where K*° and p° are completely transversely 
polarized, implying a large P-wave which is compensated by the longitudinal S'-wave. Recall 
that the factorizable D — > UU amplitudes have the sailent feature : 



IS— wave amplitude) > |P— wave amplitude| > \D— wave amplitude|. 



(45) 



Since the color-suppressed P-wave amplitude of P° — > K*°p° is proportional to [1 + ci/ (c2 + 
|ci)(A2 /A2 P )], a large P-wave thus indicates a negative A% /A2 and 



A" 



<< 



or 



A? 



0, 



7^0. 



(46) 



Therefore, we see that nonfactorizable terms in charm decay are consistently to be negative 
[cf. Eqs.(12) and (44)]. Unfortunately, at this point we cannot make a further quantitative 
analysis due to unknown final-state interactions and lU-exchange contributions. A measure- 
ment of helicities in D° — > K*°p°, D + — > <pp + will be greatly helpful to pin down the issue. 
In particular, the color- and Cabibbo-suppressed mode D + — > 0p + is very ideal for this 
purpose since it is not subject to FSI and nonspectator effects. A polarization measurement 
in this decay is thus strongly urged (though difficult) in order to test if A^ /A 2 plays a more 
essential role than A^ jA\ in charm decay. 



5. Discussion and conclusion 

The factorization assumption for hadronic weak decays of mesons can be tested on two 
different grounds: (i) to extract the effective parameters ai and especially a 2 from M — > 
PP, VP decays to see if they are process independent, and (ii) to measure helicities in 
M — > VV decay. Using the q 2 dependence of form factors suggested by QCD-sum-rule 
calculations and by some theoretical arguments, we have updated our previous work. It is 
found that a 2 is evidently not universal in charm decay. The parameter x 2 , which measures 
the nonfactorizable soft-gluon effect on the color-suppressed deacy amplitude relative to the 
factorizable one, ranges from — | to —0.60 in the decays from D — > Kn to D + — > 07r + , D — > 
K*n. By contrast, the variation of a 2 in B — > ipK, B — > D(D*)7i(p) is negligible and 
nonfactorizable terms X2(B — > ipK), Xi{B ~^ D^n(p)) are of order 10% with a positive 
sign. The pattern for the relative magnitudes of nonfactorizable effects 

\ X (B - PP, UP) I < \ X (D ^PP)\< \ X (D ^ UP) I 
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is thus well established. This means that nonperturbative soft gluon effects become more 
important when the final states are less energetic, allowing more time for final-state interac- 
tions. This explains why a 2 is class (PP or VP mode) dependent in charm decay, whereas 
it stays fairly stable in B decay. 

Taking factorization as a benchmark, we see that the nonfactorizable terms necessary 
for describing nonleptonic D and B decays are in opposite directions from the factorization 
framework. On the one hand, the leading 1/N C expansion, which amounts to a universal 
X — — §, improves the naive factorization method for charm decays. On the other hand, 
the naive factorization hypothesis works better than the large-iV c assumption for B decays 
because nonfactorizable effects are small, being of order 10%. The fact that x is positive 
makes it even more clear that the large- N c approach cannot be extrapolated from D to B 
physics. Theoretically, the next important task for us is to understand why x is negative in 
D decay, while it becomes positive in B decay. 

As for M — > VV decay, a priori effective ai )2 cannot be defined since, as pointed out 
by Kamal and Sandra, its amplitude (factorizable and nonfactorizable) involves three in- 
dependent Lorentz scalars, corresponding to S, P and D waves. This turns out to be a 
nice trade-off for solving the puzzle with the large longitudinal polarization fraction Yl/Y 
observed in B — > ipK*, which cannot be accounted for by the factorization hypothesis 
or by nonfactorizable effects weighted in the same way in all three partial waves, namely 
Ai f /Ax = Af/A 2 = V nf /V. A large Y L /Y can be achieved if B -> ipK* is almost all S- 
wave, implying that nonfactorizable contributions are dominated by the ,S-wave. The same 
assumption is also needed for understanding the ratio B(B~ — > D*°p~)/B(B° — > D* + p~). 
We found that all nonfactorizable terms A\ f /Af K *, A\ s /A^ p , A\ s /Af D * are of order 10% 
consistent with X2(B — > D(D*)n(p)) and X2(B — > ^K). 

Surprisingly, the assumption of S-wave dominated nonfactorizable effects is not operative 
in D — > K*p decay, which exhibits again another disparity between B and D physics. We 
found that A%* /A 2 should play a more pivotal role than A^ /A 1 in charm decay. We thus 
urge experimentalists to measure helicities in the color- and Cabibbo-suppressed decay mode 
D + — > (j)p + decay to gain insight in the nonfactorizable effects in D — > VV decay. 



ACKNOWLEDGMENT 

This work was supported in part by the National Science Council of ROC under Contract 
No. NSC84-2112-M-001-014. 



14 



REFERENCES 



1. See e.g. H.Y. Cheng, Int. J. Mod. Phys. A4, 495 (1989). 

2. N. Deshpande, M. Gronau, and D. Sutherland, Phys. Lett. 90B, 431 (1980); M. 
Gronau and D. Sutherland, Nucl. Phys. B183, 367 (1981). 

3. H.Y. Cheng, Phys. Lett. B335, 428 (1994). 

4. M. Wirbel, B. Stech, and M. Bauer, Z. Phys. C29, 637 (1985). 

5. A.J. Buras, J.-M. Gerard, and R. Riickl, Nucl. Phys. B268, 16 (1986). 

6. CLEO Collaboration, M.S. Alam et al, Phys. Rev. D50, 43 (1994). 

7. A.N. Kamal and A.B. Santra, Aberta Thy-31-94 (1994). 

8. M. Gourdin, A.N. Kamal, and X.Y. Pham, Phys. Rev. Lett. 73, 3355 (1994). 

9. R. Aleksan, A. Le Yaouanc, L. Oliver, O. Pene, and J.-C. Raynal, DAPNIA/SPP/94- 
24, LPTHE-Orsay 94/15 (1994). 

10. J.M. Soares, TRI-PP-94-78 (1994). 

11. M. Bauer, B. Stech, and M. Wirbel, Z. Phys. C34, 103 (1987). 

12. Q.P Xu, Phys. Lett. B306, 363 (1993). 

13. C.A. Dominguez and N. Paver, Z. Phys. C41, 217 (1988); A. A. Ovchinnikov, Sov. 
J. Nucl. Phys. 50, 519 (1989); Phys. Lett. B229, 127 (1989); V.L. Chernyak and 
I.R. Zhitnitski, Nucl. Phys. B345, 137 (1990); S. Narison, Phys. Lett. B283, 384 
(1992); V.M. Belyaev, A. Khodjamirian, and R. Riickl, Z. Phys. C60, 349 (1993); P. 
Colangelo, BARI-TH/93-152 (1993). 

14. P. Ball, V.M. Braun, and H.G. Dosch, Phys. Lett. B273, 316 (1991); Phys. Rev. D44, 
3567 (1991); P. Ball, Phys. Rev. D48, 3190 (1993). 

15. P. Colangelo and P. Santorelli, Phys. Lett. B327, 123 (1994). 

16. K.C. Yang and W-Y.P. Hwang, NUTHU-94-17 (1994). 

17. M. Witherell, in Proceedings of the XVI International Symposium on Lepton-Photon 
Interactions, Ithaca, 10-15 August 1993, eds. P. Drell and D. Rubin (AIP, New York, 
1994). 

15 



18. L.L. Chau and H.Y. Cheng, Phys. Lett. B333, 514 (1994). 

19. Particle Data Group, Phys. Rev. D50, 1173 (1994). 

20. Mark III Collaboration, J. Adler et al, Phys. Lett. B196, 107 (1987). 

21. H.Y. Cheng and B. Tseng, IP-ASTP-21-94, |hep-ph/9409lQB] (revised) (1994). 

22. I. Halperin, TECHNION-PHYS-94-16 (1994). 

23. B. Blok and M. Shifman, Nucl. Phys. B389, 534 (1993). 

24. ARGUS Collabotation, H. Albrecht et al, DESY 94-139 (1994). 

25. CDF Collaboration, FERMILAB Conf-94/127-E (1994). 

26. M. Gourdin, YY. Keum, and X.Y. Pham, PAR/LPTHE/95-01 (1995). 

27. C.E. Carlson and J. Milana, WM-94-110 (1994). 

28. Mark III Collaboration, D. Coffman et al, Phys. Rev. D45, 2196 (1992). 

29. E691 Collaboration, J.C. Anjos et al, Phys. Rev. D46, 1941 (1992). 



16 



